TO BE SUBMITTED TO ASTROPHYS. J. 

Preprint typeset using I^'T^]X style emulatcapj v. 5/2/11 



(N 



GENERIC MODEL FOR MAGNETIC EXPLOSIONS APPLIED TO SOLAR FLARES 

D. B. Melrose 

Sydney Institute for Astronomy, School of Physics, The University of Sydney, NSW 2006, Australia 

to be submitted to Astrophys. J. 

ABSTRACT 

An accepted model for magnetospheric substorms is proposed as the basis for a generic model for 
magnetic explosions, and is applied to solar flares. The model involves widely separated energy- 
release and particle-acceleration regions, with energy transported Alfvenically between them. On a 

■ global scale, these regions are coupled by a large-scale current that is set up during the explosion by 

■ redirection of pre-existing current associated with the stored magnetic energy. The explosion-related 
I current is driven by an electromotive force (EMF) due to the changing magnetic flux enclosed by this 
. current. The current path and the EMF are identified for an idealized quadrupolar model for a flare. 
I Subject headings: Sun: flares; Magnetic fields; Magnetic reconnection; Acceleration of particles 

I L INTRODUCTION 

■ Recent dramatic improvements in observational resolutions of the active Sun have led t o imp roved understa nding 
[ and modeling of detailed phenomena in solar flares ()Benz fc Giidell 120101 : iHolman et al\ 120111 : iKontar et al\ 120111 : 

• • iWhite et aLll20lil ). However, this emphasis on details has left unresolved several long-standing problems relating to 

i-C ■ the global electrodynamics and energetics of flares. We do not have an accepted global model for a flare that allows a 
semi-quantitative estimate of the power and energy in terms of large-scale observable parameters. A global model for 

Q . solar ilares is needed not only for the sola r application, b ut also as the basis for a generic model for other magnetic 

$H ■ explosions, such as flares on other stars ()Benz fc Giidell [2010l ) and magnetar outbursts ()Lvutikovll 2006[) . For any 

] form of magnetic energy release, integration of the source term (J ■ E) for electromagnetic energy over the volume in 

Ci ■ which the energy is stored implies that the power released can be identified as /$. A global description of a magnetic 
explosion requires some form of circuit model in which an electromotive force (E MF), driv es a current, /, around a 

^ ] circuit. Despite interest in circuit models in the 1970s and 1980s (Col gate 1 1978t fs"picei]ll982[ ). we have no acceptable 

■ model (a) for the circuit that describes how the stored energy is released, (b) for the current (intrinsically) associated 
^ ^ • with the energy release, nor (c) for the EMF that drives this current. Another long-standing "number problem" 

\^ , (jBrown et aZ.1 120091 ) is that the X-ray data imply energetic (e w 10^ eV) electrons precipitating at a rate N such that 
CO ' eN exceeds any plausible current by a very large factor, and such that the total number of electrons accelerated greatly 
, exceeds the number initially in the flaring coronal flux tubes. A further loi ig-standing problem is that any estimate 
. of $ based on the changing magnetic field implies a value, of order 10^° V (jSwanniri933t lColgatelll978| ). that greatly 
' exceeds e/e w 10^ V, and for which there is no direct evidence, 
psj . There is a long history of relating models for solar flares and models for magnetos pheric substorms, wh ich are a 
I ■ terrestrial form of magnetic explosion. For example, the widely favored CSHKP model (jForbes fc Malhereb6.1986 ) for 

■ flares involves a current sheet in an outflowing plasma that is analogous to the current sheet in the Earth's magnetotail. 
A major difference between flares and substorms, from a global viewpoint, is that the stored magnetic energy in a 
flare must be associated with currents that flow through the photosphere and close below it, whereas the stored 

; ^ ■ magnetic energy in a substorm is associated with currents confined to the magnetotail. From an energetic viewpoint, 
] a global mo del for a flare must involve a recon figuration of the coronal paths of currents that flow through the 
photosphere (|Melroselll997l : llTardv. Melrose fc Hud son 1998), rather than a change in the magnitu de of these currents . 
This essential ingredient is not compatible with a stand-alone CSHKP model. It has been argued (jUchida et aLl[T999l) 
that the CSHKP model should be regarded as one part of a global model for a flare, rather than as a global model on 
its own. 

In a magnetic explosion, magnetic energy builds up in a diffuse plasma over a relatively long time, and is released on a 
much shorter (explosive) time. A substantial part of the released energy appears in fast electrons (first-phase electrons 
in a flare and auroral inverted-V electrons in a substorm). In formulating a generic model for magnetic explosion, 
it is important that the model account both for the features common to flares and substorms, and for essential 
differences between flares and substor ms. There is a widely-accepted global model for substorms (jMcPherron et aU 
119731 : iTreumann &: BaumiohannI I1997D . illustrated in Fig. [T] some features of this model should to be included in 
a generic model (and applied to flares), while other features are specific to substorms. Three specific ingredients 
are assumed here to be generic. First, the primary energy release site is far removed from the region where the the 
acceleration of electrons occurs. Second, a new current system is set up during the explosion, and this current connects 
the energy-release and acceleration regions. Third, the energy transport between the energy-release and acceleration 
regions is Alfvenic, involving the parallel current and a cross-B electric field. One feature that is favored for substorms 
but cannot apply to flares concerns the driving force: this attributed to a pressure gradient in the generator region in 
Fig. [1] but pressure gradients in the corona are inadequate to drive a flare, and the driver must be an electromagnetic 
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Fig. 1. — Model for a substorm: the driving force (a pressure gradient) is indicated in the x-y plane at z = zq in the magnetotail; the 
cross-tail current is redirected along field lines, c losing in th e ionosphere (labeled Aurora) at z = zj; the dissipation occurs in a separate 
region, z = Za, between these two regions. [After iRonnmarQ ((2002 ')] 

(Maxwell) stress. Conversely, the generic model suggests that a more plausible driver for a substorm is a Maxwell 
stress analogous to the the driver for a flare. 

In the generic model proposed here, emphasis is placed on an essential aspect of the electrodynamics that is obscured 
in most models for solar flares: the intrinsic time dependence. Flares are obviously intrinsically time-dependent 
phenomena, implying an intrinsically time-dependent magnetic field. Nevertheless, most models for flares are steady- 
state: models for magnetic reconnection, electron acceleration, etc., have the intrinsic time-dependence replaced by 
proxies involving boundary conditions that allow inflow or outflow. Any steady-state model has d'Q/dt = 0, so that 
there is no inductive electric field (one with curlE ^ 0) and no displacement current {dE/dt = 0); the electric field in 
the model is necessarily an electrostatic or potential field. The physical differences between an inductive field and a 
potential field are ignored in most models. For some purposes, the steady-state models are effective proxies. However, 
this is not the case for other purposes, notably for electron acceleration. Although a variety of alternative acceleration 
mechanisms remain under consideration, the see mingly simplest and most plausibl e mechanism, acceleration by a 
parallel electric field, , is so poorly understood ()Melrosell2009l : iBenz fc Gudel|[20Tot ) that it is not always considered 
as a viable mechanism for flares (jZharkova et a/.ir201lD . whereas it is the accepted mechanism for the acceleration of 
aur oral electrons in a su bstorm. Inclusion of explicit time-dependence is essential for understanding acceleration by 
£^11 (|Song fc LvsakI 120061) . 

The generic model for magnetic explosions proposed here has one additional notable feature. Besides being based 
on substorms and including time-dependence explicitly, it involves an explicit separation into multiple scales. Three 
scales are introduced: global, macro and micro. The global scale is at least as large as an active region, and is the scale 
of the circuit around which the flare-associated current flows. The macro scale includes the energy-release region and 
the acceleration region, and the energy transport between them. The collisionless plasma processes that are essential 
for magnetic reconnection (in the energy-release region) and for E'n (in the acceleration region) occur on the micro 
scale, related to the Debye length and other fundamental lengths in the plasma. Such collisionless processes are highly 
localized and transient, and a statistically large number of them is required to have macroscopic consequences. In 
this paper, emphasis is placed on the global scale. Detailed models for some of the processes on the macro scale are 
discussed in an accompanying paper (Melrose 2012), referred to as paper 2. 

The generic model for a magnetic explosion is summarized in Sj2] Application of the model to solar flares, involving 
identification of the current and the current path, is discussed in !j3l and the EMF is discussed in §4l Comparison with 
the application to substorms is discussed in Sj5l The results are discussed in 5^6] 

2. GENERIC MODEL FOR MAGNETIC EXPLOSIONS 

In this section a widely accepted model for magnetospheric substorms is reformulated as a generic model for magnetic 
explosions in space and astrophysical plasmas. 

2.1. Generic model 

The generic model combines several ideas: nonlocal energy release, a multi-scale approach, intrinsic time-dependence, 
and globally controlled dissipation. The generic model contains the following ingredients: 

1). Driver:: The temporally changing magnetic field produces an EMF, $. 
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2) . Energy-release region:: The EMF drives a current, /, across magnetic field lines in an energy-release region and 

around a (new) global-scale circuit. 

3) . Energy transport:: The magnetic energy released is transported away from the energy-release region as an Alfvenic 

flux. 



4) . Current closure:: The (newly directed) current flows along field lines to a remote closure region, and back to the 

energy-release region along neighboring field lines. 

5) . Acceleration region:: The current fiows through an acceleration region, where $ partly localizes along field lines, 

leading to electron acceleration by and dissipation of the released magnetic energy. 

6) . Globally controlled dissipation:: CoUisionless dissipation processes on a micro scale combine to give an effective 

dissipation whose rate is governed by the global requirements. 



2.2. Remarks on the generic model 
Several remarks on these six ingredients in the generic model are appropriate. 

1). The words "driver" and "trigger" are used with specific meanings in the context of magnetic explosions. There 
are two different drivers. One applies to the build-up phase, is different for substorms and flares, and is not directly 
relevant to the explosion. The effect of this driver is to store magnetic energy in a metastable conflguration on a 
relatively long time scale. The trigger for a magnetic explosion initiates the (rapid) change in configuration that 
results in release of the stored energy; the trigger is not clearly identified in any context. The primary magnetic energy 
release requires magnetic reconnection, and the "driver " referred to here drives infl ow into a reconnection region. The 
only possibilities for this driver are a pressure gradient (|Satolll978l : lHaerendeill2007f ) and the Maxwell stress (the J x B 
force in MHD). The assumption that a pressure gradient is the driver is questionable for a substorm, and unacceptable 
for a flare. When the intrinsic time dependence is taken into account, the Maxwell stress becomes a natural driver in 
both cases, as discussed in detail in paper 2. 

The name "generator region" is used in the magnetospheric literature for the energy-release region, and this name 
is regarded as misleading in the present context. In both substorms and flares, pre-existing currents are redirected, 
rather than an intrinsically new current being generated. 

The energy release involves magnetic reconnection. However, an essential aspect that is not addressed in most 
detailed models for solar flares is that the reconnecting fleld is twisted, and carries a current that is appro ximately 
force - free. It is the redirection of this current on the global scal e that releases the stored magnetic energy ([Melros' 



1997f). Some detailed calculations f or reconnection are available (jPahlburg. Antiochos fc Nortonlll997L iLintonf 



200 



Longcope. Guidoni &: Lintonll 20091 ). Here it is simply assumed that magnetic flux and current are transferred together 



during reconnection. 

2) . Energy flow is determined by the Poynting vector, with E — — u x B in ideal MHD: 

S^^=\^u^B. (1) 

A*o Aio Mo 

The inflow of energy into the energy- release region may be described by the flrst term on the right hand side of ([1}, 
with u-B — 0; the outflow may be described by the second term, referred to here as an Alfvenic flux. The argument for 
Alfvenic energy transpo rt has b een made in connection with s everal flare models (jEmslie fc Sturroc3ll982t iHaerendel 
119941 : [Fletcher fc Hudso n 2008; Kontar. Hannah fc Bianir201lD . In other flare models the energy outflow is assumed to 
consist of energetic particles and/or mass motions. 

3) . Alfvenic energy transport involves a parallel current, a cross-B electric fleld, and a vortex motion, which can 
be interpreted as a propagating magnetic twist. Process es whereby the magnetic en ergy released by reconnection in a 
flare is converted into an Alfvenic flux were discussed by [Fletcher fc Hudson! (|2008f ) . An idealized cylindrical model is 
presented in paper 2. 

4) . In a substorm, parallel currents flow up and down along flux tubes, which converge just above the ionosphere, 
where the current closes across field lines. This closure is usually attributed to the Pedersen conductivity, which can 
be described by an effective resistance, Rp, determined by the inverse of the height-integrated conducti vity. In a flare , 
the current may close in a similar manner in a partially ionized region of the chromosphere (Melrose fc Khanlll989f ). 
Alternatively, the current may close dynamically, driving an Alfven wave into the denser regions of the atmosphere 
where the Alfvenic propagation time is long compare d with the time scale of the explosion, due to the rapid increase 
of density with depth through the chromosphere (Wh eatland fc Melrosel[T994[) . The effective resistance of the closure 
region can be interpreted as Rp or the Alfvenic impedance, Ra = fJ-o'^^A'- Alfvenic energy propagation dominates for 
Ra ^ Rp, and resistive dissipation of the energy dominates for Rp ^ Ra- In the model assumed here, the impedance 
of the current-closure region is smaller than other relevant impedances, and the energy dissipation associated with this 
current closure is negligible. 

5) . Electrons are assumed to be accelerated by £'||, which develops by $ becoming partly localized along the field 
lines. This acceleration process is inadequately understood, as discussed in paper 2. The energy released goes into bulk 
energization of electrons in compact flares, which are non-eruptive. In an eruptive flare, much of the energy released 
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goes into kinetic energy of mass motion of a coronal mass ejection (CME). The acceleration of the CME requires a 
modification of the generic model that is not discussed here. 

6). CoUisionless dissipation involves a statistically large number of localized and transient instabilities. This allows 
a major simplification: one can assume that the statistical distribution adjusts to meet the global requirement that 
drives the dissipation. In a driven system, the dissipation is maximized when the dissipation corresponds to an effective 
resistance, i?cff of order Ra, when the energy is absorbed, rather than being transmitted (for i?eff ^ Ra) or reflected 
(for Rcs ^ Ra)- The effective dissipation associated with the conversion to an Alfvenic flux in the energy- release 
region and with the actual dissipation in the acceleration region are assumed to correspond to RcS ~ Ra- 

3. FLARE-ASSOCIATED CURRENT 

In the application of the generic model to solar flares, unlike that to substorms (^S)), the essential features of the 
model do not follow directly from observations. On the global scale, the essential features that need to be identified 
are the current and the EMF. The current is discussed in this section and the EMF in fjH 

3.1. Release of magnetic free energy 

Prior to a flare, magnetic ene rgy is stored in a complex magnetic structure overlying an active region. Newly 
emerging flux tubes are twisted (jLeka et aZ.1 119960 . implying that they carry currents. As the magnetic flux builds 
up, the current builds up correspondingly. A tw isted flux loop has a current flowing from one footpoint to the other 
and closing at the based of the convection zone (jMcClvmont fc Fis her 1989). Magnetic energy in the corona may be 
attributed to three different classes of current: (a) currents that are confined to the solar interior, (b) currents that 
are confined to the corona, and (c) currents that link the interior and the corona. Currents of class (a) produce a 
"potential" component of the coronal magnetic field that is unchanged during a flare. The energy stored in currents 
of class (b) is too small to account for the energy released during a flare; specifically, the assumed closure in the 
corona implies J x B in the closure region, and this force must be balanced by a mechanical force (a pressure gradient 
or gravity), which severely restricts the magnitude of J x B. Although these currents may be important during a 
flare, when the J x B be balanced the acceleration of plasma flows, they cannot be important energetically on the 
global scale. On ly currents of class (c), wh ich flow along fleld lines in the corona, and close near the base of the 
convection zone (IMcClvmont fc Fished Il989l) are relevant to the global energy storage. The current passing through 
the photosphere/chromosphere can be determined from vector magnetograms, which imply currents in flaring flux 
tubes > 10^^ A, and under extreme conditions « 10^'^ A (jMoreton fc Severnvlll968l ). 

The magnetic free energy is associated with the resulting current system, with all (energetically) important currents 
flowing through the photosphere. On a global scale, the free energy can be identified as \LI^, where L = ji^l is an 
effectiv e inductance, and i h as the dimensions of a length. The net current cannot change at the footpoints during 
a flare (jMelrosd 119951 Il997| ). implying that I does not change during a flare. Except for regions around magnetic 
nulls, where reconnection occurs, these large-scale coronal currents flow along fleld lines, and a change in the current 
conflguration must be associated with a change in the magnetic conflguration. 

The trigger for a flare has not been identifled explicitly. The newly emerging flux builds up magnetic stresses in the 
corona over a time scale of days, and the reconnection must allow a reconfiguration of the magnetic field and current 
on a large scale to account for the magnetic energy release. In an emerging-flu x model, the trigge r presumably involves 
the emerging flux tube being forced into contact with an overlying flux tube (jNishio et a/.lll99"7l ). 

In summary, the release of magnetic free energy in a flare requires that magnetic reconnection in the corona leads 
to a magnetic reconflguration that reduces the effective length, I = L/ ^q, of the current-carrying flux tubes. 

3.2. Flare models 

Although widely favored, the CSHKP model (Fig. [2]) is inadequate as a global model for a flare. In order to discuss 
the global energetics, one needs a model that relates the currents in the corona to the currents passing through the 
photosphere, and the global current is not included in a conventional CSHKP model. This deficiency in the CSHKP 
model was one of the motivations for developing a quadrupolar model. 

A quadrupolar model is the simplest example of a magnetic (and current) configuration t hat allows ma gnetic en- 
ergy release without changing the current or magnetic fiu x at the photosphcric boundary (Melrose" 1997). Various 
qua drupola r mod els have been proposed (Uchida 1996; Nis hio et all 1997; Uchida et al. 1999; Aschwandcn et al. 199!|; 
lAchwandenl 1200^ . One version is an idealized form of an emerging fiux model, in which a rising magnetic loop 
encounters an overlying magnetic loop. As in Fig. |3l let the two loops be labeled 1 and 2 with footpoints 1± and 
2±, respectively. Reconnection is assumed to occur at the point of contact, C say, forming two new flux loops con- 
necting l-i- to 2_ and 2-i- t o 1_ respectively. Numerical models for reconnection of tw o current-carrying flux loops 
(|Linton fc Longcopil 20061 ) have recently been applied to a quadrupolar configuration (jTorok et a?.ll20HT ). showing 
how the two new fl ux tubes are formed. An analytic model allows the magnetic energy released to be estimated 
(jHardv. Melrose fc" H udson 1998) in terms of the magnetic reconfiguration reducing the effective length of the current 
path. Such models describe the initial and final states, with the energy released being identified as the difference in 
the stored magnetic energy between these two states. The current during the flare, indicated schematically in Fig. ^p, 
is discussed separately below. 

A more realistic conflguration for many flares involves an arcade of current-carrying magnetic loops. A shearing 
motion that displaces the footpoints in opposite directions on either side of the arcade increases the length of the 
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Fig. 2. — The standard CSHKP model for a flare as illustrated bv lForbes fc Malherebd l l 19861 ). In the interpretation assumed here, the 
reconnecting flux tube carries a current, and the energy release is due to the shortening of the current path as a reconnected flux tube 
moves closer to the photosphere. 

current paths, storing magnetic energy. Magnetic reconnection between neighboring flux tubes, which pairwisc form 
a quadrupolar-type configuration, can lead to a net shortening of the current path, and release of the stored energy. 

The shortening of the current path has potentially observable consequences. While the difficulties in identifying the 
change in magnetic configuration from pre-flare to post-flare are well know n, vector m agnetogram data suggest an 
increase in the transverse component of the photospheric field during a fiare (|Wang|l2011f ). This is consistent with the 
shortening of the current path: qualitatively, a shortening of the net path requires a reduction in the angle at which 
the magnetic field emerges through the photosphere, implying an increase in the tangential component of the field. 



The new current path set up during a flare must have the effect of converting the initial current configuration into 
the final current configuration. A current that leads to this result is found by subtracting the initial current from the 
final current. 

For the quadrupolar model illustrated in Fig. [3b, the difference between the final and initial current paths separates 
into two paths. Let the point where the two initial current paths (thick solid lines) reconnect at the point of contact, 
C, be identified as the energy-release region. These two paths (projected onto the solar surface) are the triangles 
2_ C1+ and 2+ 1_ C2_|_. Let the initial currents between 1— and 2-|-,2— be 1\ and /2, respectively, and the 
final currents between these footpoints be J{ and respectively. Conservation of current at each footpoint requires 
I[ — Ii — II and I2 = I2 — II, with II the final current between l-f,2— and 2+, 1 — . The fiare-associated current 
in both triangles is I^. The current closes across field lines at each of the four footpoints and at C. Thus, there 
are up and down fiare-associated currents at all the four footpoints where the initial and final loops meet at the 
photosphere/chromosphere. Downward electron acceleration, producing hard X-rays, is expected with the upward 
current at each of the four footpoints, by analogy with the auroral regions in substorms. 



Although the quadrupo lar model invo lves global parameters / and it is not a circuit model in the conventional 
sense. In a circuit model ()Colgatelll978l : FSpicer 1983), the current path is specified, and the EMF and the power are 
associated with a photospheric dynamo. In the model illustrated in Fig. [3l the currents flow along field lines, except 
at the point C and the closure regions at the footpoints. During the flare, the current configuration changes due to 
transfer of current-carrying fiux tubes from the initial to the final flux loops; the light solid line in Fig. [3}d indicates an 
intermediate location of a flux tube in motion. The intrinsic time dependence of the model is an essential ingredient is 
estimating the EMF $ (2]). There is no photospheric dynamo. Although the physics is quite different, the numerical 
values of / and $ are similar in the two models. These parameters are constrained by the requirement that power, /<&, 
be equal to that released in a flare, and that the current be of order that inferred from magnetograms, / = 10^^-10^'^ A. 



3.3. Current path during a flare 



3.4. Comparison with circuit models 
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2+ 1+ 2+ 

Fig. 3. — Idealized quadrupolar model for a solar flare: (a) initial loops in gray, flnal loops in black; (b) projection onto the solar surface 
with the initial loops shown as thick solid lines, the final loops as dashed lines, and reconnected flux loops moving from one to the other, 
shown by thin lines. The EMF is assumed to dr ive current around the two triangles in (b), requiring current closure across field lines at 
each of the four footpoints. [After IMelrosel {1993)] 

A more subtle difference concerns the dissipation. In circuit models the energy release it attributed to Ohmic 
dissipation, described by a resistance in the circuit. The quadrupolar model is intrinsically time dependent, and the 
conventional circuit equation for energy does not apply to a time-dependent model. One can use energy arguments 
only by comparing the initial and final states, which are assumed time-independent. The EMF, $, is (minus) the rate 
of change of the magnetic flux, and the counterpart of a resistance is played by (minus) the time derivative of the 
effective inductance, which decreases as the effective current path length decreases. 

4. FLARE-ASSOCIATED EMF 

The EMF that drives a flare must be associated with the changing magnetic field. In this section, three idealized 
models for a time- varying magnetic field are used to estimate the EMF. 

4.1. EMF associated with magnetic energy release 

The first estimate of $ is based on the integrated form of curlE = —d'B/dt applied to the time- varying magnetic 
field within the energy-release region. 

Consider a model in which the time- varying magnetic field is confined to a box with sides of length Lx,Ly,Lz, that 
is, the region < x < —Ly/2 < y < Ly/2, —Lz/2 < z < Lz/2. Let the magnetic field be along the z axis 

with 

B,{x,t) = Boi2x/L,)f{t), (2) 

where f{t) is a decreasing function over a characteristic time T, e.g., f{t) = exp(— i/T) or 1 — t/T. The inductive 
electric field is along the y axis, with 

Ey{x, t) - Ey{0, t) = f{t)Bo{xyL,), (3) 

where the dot denotes the time derivative. The potential, ^{x,t) = ~f{t)Bo{x'^Ly/Lx), across the box along the 
y axis, is symmetric on either sides of the plane x = 0. The current density (obtained by integrating the induction 
equation ignoring the displacement current) is along the j/-axis with Jy{t) = —Bo{2/iJ,oLx)f{t)^ implying a current 
I{t) = —2BoLzf{t)/iiQ independent of x. The magnetic flux has one sign for a; > and the opposite sign for x < 0, 
with 'i>B — iBoLxLyf{t)/2 in the two halves. The rate of change of the magnetic energy, dWB/dt, is given either 
by the time derivative of Wb, obtained by integrating the magnetic energy density, B^{x,t) /2^o, over the box, or by 
averaging I{t)^{x,t) over the box, either of which gives dWB/dt = f(t)f(t)BQLxLyLz/&Ho- 

This model implies an EMF of order $ — BoLxLy/T, where T is the duration of the fiare. Fiducial numbers are 
discussed below, and one choice is Bq — 10^^ T, — Ly — 10^ m, T = 10^ s, imply ing $ = 10^" V. That such a large 
potential is needed to drive a fiare has been known for many decades ()Swan 

The neglect of the displacement current in the foregoing calculation is well justified for [Lx/cTY ^ 1- However, 
as discussed in paper 2, the displacement current drives a polarization current that is a factor /v\ larger, and the 
condition for the neglect of this current, [Lx/vaTY ^ 1, is not well justified. It is argued in paper 2 that the cross-field 
current associated with the energy release is the polarization current, and hence that the displacement current plays 
an essential (albeit indirect) role in the energy release. 

4.2. Time-varying inductance 

The second estimate of the EMF is based on the rate of change of the magnetic flux enclosed by the circuit around 
which the current flows. As argued above, the change in stored magnetic energy can be attribute to a reduction in 
the effective inductance, L{t) — /ioi!(i), of the circuit, requiring a reduction in the effective length, €(t), of the current 
path. The magnetic flux in such a model is 

^B{t)^L{t)I, (4) 
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implying that the EMF is 

= -L{t)I = -fioiit)!. (5) 

One can write this in the form of Ohm's law, — J^o{t)I, with an effective impedance So(0 = fJ-()Vo{t), vo{t) — ~£{t). 
Further evaluation of this model requires an estimate of £, which depends on the geometry of the circuit. 

Changes in the stored magnetic energy due to changes in the geometry are ir nportant in laboratory plas ma ex- 
periments, and can be modeled through changes in self and mutual inductances ([Book. Turchi fc SteinI 119781 ) . Any 
distribution of currents in the corona can be approximated by a set of sub-currents, I = Ii + I2 + • ■ • say, such that 
the magnetic energy can be written as the sum i(Li/^ + L2I2 +■■•) + -^^12-^1 ^2 + • • • : where Li,L2, ■ ■ ■ are self 
inductances and M12, • • • are mutual inductances. This sum is written as ^LI^ to define L. When the sub-currents are 
close together, the self and mutual inductances are nearly equal (to L{t)) and the mutual inductances decrease with 
increasing separation between the sub-currents. Ass uming that a single value of L{t) suffices, simple models for the 
geometry ()Book. Turchi &: SteinllTOTSt lMelroselll997fl imply that £ differs from the actual length of the current path by 
a factor of order unity, which can be ignored for semi-quantitative purposes. One can then estimate vq — A£/T, where 
A£ is the difference between the initial and final lengths of the current path. Thus this estimate gives $ = fioA£I/T. 
The fiducial numbers discussed below suggest A£ — 10^ m, / = 10^^ A, T — 10^ s, implying $ ~ 10^° V, consistent 
with the first estimate. 

4.3. Propagating flux tube 

The third estimate of $ is based on transfer of magnetic flux and associated current from the initial to the final 
flux tubes. Reconnection of force-free flux tubes in 3D leads to a portion of the initial flux tube propagating from 
its initial location to its final location, as illustrated in the quadrupolar model. Fig. The effective length, £{t), 
of the flux tube reduces, at the rate vo(t), as it propagates. This model gives an alternative way of estimating 

^{t) = —L{t)I — —fiQVo{t)I by providing a geometric model for vo{t). 

Consider either of the triangles in Fig.[3}D, which are projections onto the solar surface. This suggests approximating 
the reconnecting flux tube as a triangle. Assuming that the final flux tube is in a vertical plane, whose base is the 
dashed line, the initial flux tube has the same base, and is in azimuthal plane, at an angle A(j) say to the vertical. Let 
the base be of length £0, and the initial and final triangles being equilateral with sides ^1/2,^2/2, respectively, with 
perpendicular bisectors of length di,d2, respectively. A flux tube in the process of being transferred is at < (?i < A0 
with sides £/2 and perpendicular bisector d. It is convenient to introduce the angle C, such that £ = 2ds'm(. One has 
£ = 2dsi'D.(. The velocity of the apex of the triangle may be identified as Vq; it has contributions from the rate of 
change of both ( and (j), although only the former contributes to £. To within a factor of order unity, one may estimate 
Vq d ^ £. With £ = ~A£/T, this model reproduces the results of the other two models. 

4.4. Fiducial numbers for a flare 

As fiducial values, consider a moderately powerful flare with an energy /$ = 10^^ W transferred to energetic electrons 
over a time T ~ 10^ s, due to a current 10^^ A and an EMF $ — 10^^ V. These numbers are consistent with a magnetic 
field B = 10^^ T in a box with sides £ = 10^ m, with the reduction in the effective length, A£, of the current 
path being of order 10^ m. Assuming the implied effective impedance, 10~^ J7, to be the Alfvenic impedance implies 
VA = ms~^ in the energy-release region. These numbers are also consistent with the stored energy being identified 
as ■^LI^ = 10^'^ J, with an inductance L = ^q£ = 10 H. The actual numbers for a given flare need to be greater or 
smaller than these values to account for less and more powerful flares. 

It should be emphasized that the high value of $, of order 10"'^° V, is essential, but that there is no direct evidence 
for it. This contrasts with a substorm, where the maximum energy of auroral electrons is plausibly determined by 
e*!". In flares electrons are accelerated to 10-20 keV, which is e^/AI with M w 10^. The factor is referred to here as 
the multiplicity, M, in the sense that the rate, N, electrons are accelerated is Ml/e. The multiplicity quantifies the 
long-standing "number problem." 

5. APPLICATION TO SUBSTORMS 

Application of the generic model to substorms (|Akasofulll964l : IKivelson fc Russelilll995[ ) is discussed briefly here for 
background and comparison purposes. 

5.1. Description of a substorm 

The magnetic energy released in a substorm is stored in the Earth's magnetotail. Prior to a substorm, magnetic 
flux builds up in the magnetotail due to reconnection on the sunward side of the magnetosphere at a time when the 
direction of the magnetic field in the solar wind is opposite to the E arth's magnetic field there. Reconnection occurs 
at a rate determined by a coupling function (Gonzalez 199(1 iLopez et al. 2004). The solar wind drags the reconnected 
flux tubes around to the anti-sunward side, forming the magnetotail. In a steady state, there is slow reconnection in the 
magnetotail, with the magnetic flux being returned to the sunward side through the inner part of the magnetosphere. 
The magnetic energy stored in the magnetotail is associated with a cross-tail current that flows along the dawn-dusk 
line in the current sheet that separates the oppositely directed field lines above and below the sheet. This current 
closes (around nearly semi-circular paths) on the surface of the magnetotail above and below the current sheet. When 
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the orientation of the magnetic field in the solar wind reverses (a "southward turning"), the supply of magnetic flux to 
the tail is shut off, and the steady-state balance of stresses in the magnetotail can no longer be sustained. This is one 
plausible cause for the onset of explosive energy release in a sub storm. Th e magnetic energy release is associated with 
magnetic reconnection or current disruption in the magnetotail (|Luill20l"l[ ) . The tension in the newly formed flux tube 
accelerates plasma both toward and away fro m the Earth. The observed earthward high speed flow correlates with the 
onset of a substorm (Baum iohann. Paschmann fc: Liihi1ll99Ct |Pu et a/.l[2bl0. ). Such flow is a possible magnetospheric 
counterpart of mass flows, including jets and CMEs, associated with flares. 

During the substorm, magnetic reconnection leads to partial "dipolarization" of the magnetic field, in the sense that 
some open field lines become closed, reducing the magnetic flux in the open held lines. The cross-tail current is partly 
redirected into a field-aligned current flowing out of the energy-release region, forming the so-called current wedge, 
cf. Fig. [TJ The up and down currents from opposite ends of the cross-tail current flow along converging flux tubes, 
closing across field lines in the ionosphere. There are actually two current wedges, one in each hemisphere, so that 
there are up and down currents above the current-closure regions in both hemispheres. The magnetic energy released 
is transported Alfvenically along field lines to the acceleration region, which is located far above the ionosphere in 
the upward current region. An analogous Alfvenic energy flow is assumed to occur in flares, but the location of the 
acceleration region in a flare is not as well determined. 

In the acceleration region, energy is transferred to auroral electrons by a parallel electric held, which is highly 
structured, involving a large number of localized regions (over several Debye lengths), propagating regions with ^ 
(HE gun "^r^ll998at IM ozer fc Kletzind[l998| ). modeled as phase-space holes (jErgun a/.lll998bf) or as kinetic A lfven 
waves (jChaston et aLll2003[ ). The acceleration can be so efflcient that charge starvation occurs (|Ergun et a/]|2004|) : the 
current-carrying field lines can become devoid of thermal plasma, forming an auroral cavity. This tends to produce a 
flux of precipitating ( "inverted- V" ) electrons with energy « e$ and precipitation rate dN/dt « I /e. There are nearby 
associated downward current regions, in both hemispheres, where ions are accelerated downward and the electron 
acceleration is associated with intense electrostatic turbulence (jErgun "^r^ll998aLl2?)nl . A major difference between 
flares and substorns is that the electrons in a flare are accelerated only to a small fraction of e$. Why this is the case 
is a long-standing problem. A possible explanation is suggested in paper 2. 

5.2. Fiducial numbers for a substorm 

Substorms have an onset phase of a few minutes and a total duration of a few hours, implying T in the range 
10^-10^ T — s is adopted here. The parallel (Birkeland) currents are about 10^ A in quiet times, increasing to 
> 10^ A during a substorm. The EMF must be greater that the maximum energy of the auroral electrons. In inverted- 
V events, the energy spectrum has a maximum of several keV. This suggests flducial values J = 3 x 10^ A, $ = 10^ V. 
The magnetic flux enclosed by the circuit around which the current flows is estimated as A^^a/ = = 10^ Wb. (The 
portions of the circuit around which the Birkeland currents flow do not contribute significantly to the changing m-) 
The estimated number is consistent with a (near) magnetotail magnetic field of about 1 nT times that cross-sectional 
area of the magnetotail. 

Assuming that the internal impedance in the energy-release region is determined by Ra there, the fiducial number 
Ra — ^/I = 0.3 implies an Alfven speed, = 3 x 10^ ms~^. This is a relatively low value, which is possible for a 
region that is overdense (in a current sheet) with a weak magnetic field (around a magnetic null) . 

6. DISCUSSION 

The generic model for magnetic explosions presented in ij2]is based on a widely accepted model for magnetospheric 
substorms, and is adapted in ^and §l]to apply to solar flares. The model has separate magnetic-energy-release and 
electron-acceleration regions, with Alfvenic energy transport between them. Such a model involves different physics 
on vastly different scales, and three scales are introduced: global, macro and micro. In this paper the emphasis is on 
the global scale, which necessarily involves integrated quantities and circuit-type concepts. An obvious advantage of 
a circuit model is that it allows semi-quantitative discussion of the energy storage and release in terms of currents, as 
opposed to appealing to release of an ill-defined free magnetic energy density. However, most existing circuit models do 
not reflect the essential physics assumed here: explosive release of magnetic energy stored in the corona. Although the 
model proposed is not a circuit model in the conventional sense: the current path is not fixed, there is no photospheric 
dynamo, and the energy release is not due to the resistive decay of a current. Rather, the magnetic field is changing 
as a function of time, and this causes an EMF that redirects pre-existing current around a new path. Energy release 
is attributed to a reduction in the stored magnetic energy through a reduction in the effective length of the current 
path. Reconnection is an essential ingredient in allowing the magnetic configuration to change, to one of lower energy, 
but it plays only a secondary role in the model. 

One motivation for introducing a generic model for magnetic explosions is to avoid what I regard as a misleading 
feature in the existing literature: the use of steady-state proxies to describe intrinsically time-dependent electrodynamic 
phenomena. The inductive electric field (curl E ^ 0) and the displacement current are excluded by the steady-state 
assumption. When the time dependence is taken into account, these are nonzero and play important roles, which are 
not treated properly by the steady-state proxies. In this paper, this is shown to be the case for the EMF that drives 
the flare-associated current in a flare. Other roles that the inductive electric field and the displacement current play 
in flares are discussed in paper 2. 

One purpose of a global model for a flare is to describe the energetics. The power released can necessarily be written 
as in terms of the current, /, and the EMF, $. It is argued that appropriate fiducial values, for a flare with a 
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power 10^^ W are / = 10^^ A and $ — 10^" V. The ratio $/T, which is attributed to a resistance in a conventional 
circuit model, is interpreted as an effective impedance fJ-oVo, with vq = —i due to the decreasing length of the current 
path. It is suggested that this impedance is of order the Alfven impedance, Ra — fJ'O'^^A- Fiducial number for a flare 
suggest $// is of order 0.1 fi, which may be interpreted as a relatively slow Alfven speed, va = 10^ ms^^, in the 
energy- release region, or perhaps as vq being significantly smaller than va due to geometric factors. 
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